INFINITE DIVISIBILITY OF SOLUTIONS TO SOME SELF- SIMILAR 
INTEGRO-DIFFERENTIAL EQUATIONS AND EXPONENTIAL 
FUNCTION ALS OF LEVY PROCESSES 



P. PATIE 



Abstract. We first characterize the increasing eigenfunctions associated to the following family 
of integro-differential operators, for any a,x > 0,7 > and / a smooth function on 

(0.1) LW/(i) = x- a (^x 2 f"(x) + (aj + b)xf'(x) 

+ J™ ((f(e- r x) - f(x))e-^' + xf'(x)rl {r < 1} ) v(dr)j 

where the coefficients b G SK, a > and the measure v, which satisfies the integrability condition 
/ °°(1 A r 2 ) v(dr) < +00, are uniquely determined by the distribution of a spectrally negative 
infinitely divisible random variable, with characteristic exponent ip. L' 7 ^ is known to be the 
infinitesimal generator of a positive a-self-similar Feller process, which has been introduced by 
Lamperti [23]. The eigenfunctions are expressed in terms of a new family of power series which 
includes, for instance, the modified Bessel functions of the first kind and some generalizations of 
the Mittag-Leffler function. Then, we show that some specific combinations of these functions 
are Laplace transforms of self-decomposable or infinitely divisible distributions concentrated on 
the positive line with respect to the main argument, and, more surprisingly, with respect to 
the parameter In particular, this generalizes a result of Hartman [16] for the increasing 

solution of the Bessel differential equation. Finally, we compute, for some cases, the associated 
decreasing eigenfunctions and derive the Laplace transform of the exponential functionals of 
some spectrally negative Levy processes with a negative first moment. 



1. Introduction 

During the last decade, there has been a renewed interest for self-similar semigroups, some- 
thing which seems to be attributed to their connections to several fields of mathematics and 
more generally to many area of sciences. For instance, in probability theory, these semigroups 
arise in the study of important processes such as self-similar processes, branching processes and 
also in the investigation of self-similar fragmentation. Moreover, the Feller processes associated 
to self-similar semigroups are closely related, via the Lamperti's mapping, to the exponential 
functionals of Levy processes which appear to be key objects in a variety of settings (random 
processes in random environment, mathematical finance, astrophysics . . .). We refer to Bertoin 
and Yor [6] for an interesting recent survey on this topic. Finally, we emphasize that they are 
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also related with the theory of fractional operator which is used intensively in many applied 
fields, see e.g. the survey paper of Kilbas and Trujillo [2~T] . 

In this paper, we provide, in terms of power series, the increasing eigenfunction, associated to 
the linear operator i/ 7 ', given by (jQ.lj) . that is the increasing solution to the integro-differential 
equation, for x, q > 0, 

L (7) / 9 (^) = Qf q {x). 

As a byproduct, we compute the Laplace transform of the first passage times above for spectrally 
negative self-similar processes and some related quantities. 

Moreover, when the spectrally negative random variable has a negative first moment, we provide, 
under an additional technical condition, the decreasing eigenfunctions associated to L = 
We deduce an explicit expression of the Laplace transform of the exponential functional of some 
spectrally negative Levy processes with negative mean. This is a companion result of Bertoin 
and Yor [5] who characterized, in terms of its negative entire moments, the law of the exponential 
functional of spectrally positive Levy processes which drift to — oo. 

Furthermore, it is plain that i/ 7 ' is a generalization of the infinitesimal generator of the (re- 
scaled) Bessel process {v = and a = 2). In this specific case, Hartman |16j . relying on purely 
analytical arguments, showed, that the function 

I /^(a)Io(A) 

(1.1) 7^ T ^i , 0<a<^<oo, 

is the Laplace transform of an infinitely divisible distribution concentrated on the positive line, 
where \ v stands for the modified Bessel function of the first kind. We mention that in the limit 
case A — > oo, the result above has been reproved, in an elegant fashion, by Pitman and Yor [30] . 
We shall provide a simple probabilistic explanation of Hartman's result (for any < A < oo) and 
show that this property still hold for similar ratios of the increasing eigenfunctions associated 
to the non-local operator L^. 

The outline of the remaining of the paper is as follows. In the sequel, we set up the notation 
and provide some basic results. Section 2 is devoted to the statement of the main results. 
The proofs are given in Section 3. Finally, in the last section, we illustrate our approach by 
investigating some known and new examples. 



1.1. Notation and preliminaries. 



1.1.1. Some important set of probability measures. Let £i be a spectrally negative infinitely 
divisible random variable. It is well known that its characteristic exponent, tp, admits the 
following Levy-Khintchine representation 

POO 

(1.2) i){u) = bu + -u 2 + / (eT ur - 1 + url {r < 1} )v(dr),u > 0, 

2 Jo 

where the coefficients b S fH, a > and the measure v, image of v by the mapping x — ► —x, 
which satisfies the integrability condition (1 A r 2 ) v{dr) < +oo, are uniquely determined by 
the distribution of £i. We exclude the case when b < and J_ (1 A r) i>{dr) < +oo, i.e. when tp 
is the Laplace exponent of the negative of a subordinator. It is plain that lim^^oo ip(u) = +oo 
and by monotone convergence, one gets E[£i] = b— L rv{dr) E [— oo, oo). Differentiating again, 
one observes that V is strictly convex unless £ is degenerate, which we exclude. Note that is 
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always a root of the equation tp(u) = 0. However, in the case E[£i] < 0, this equation admits 
another positive root, which we denote by 8. This yields the so-called Cramer condition 

R[e * 1 ] = 1. 

Then, for any E[£x] £ [— oo,oo), the function u i— > ip(u) is continuous and increasing on 
[max(0, 0), oo) and thus it has a well-defined inverse function (j) : [0, oo) — >• [max(0, 0), oo) which 
is also continuous and increasing. We denote the totality of all functions ip of the form (|1.2p by 
£/C. Note, from the stability of infinitely divisible distributions under convolution and convolu- 
tion powers to positive real numbers, that CKL forms a convex cone in the space of real valued 
functions defined on [0, oo). 

We also mention that when a probability measure dm is supported on a subset of [0, oo), 
then dm is infinitely divisible if and only if its Laplace transform satisfies the conditions, for any 
u> 0, 

oo 

e~ ux dm(x) = e~^ u) 



o 



with (f)(0) = and (f>'(u) is completely monotonic, i.e. <j)' is infinitely differentiable on (0,oo) and 
for all n = 1,2..., (— l) n ~ 1 (f)( n \u) > 0, u > 0. Moreover, the so-called Laplace exponent, <f> 
admits the following Levy-Khintchine representation 



oo 



(1.3) <p{u) = au+ (1 - e' ur )fi(dr), u > 0, 

J o 

for some a > and some positive measure fj, on (0, oo) satisfying / °°(1 A r)fi(dr) < oo, see 
e.g. Meyer [27] . 

Finally, we recall that a random variable H is self-decomposable (or of class L) if it is solution 
to the random affine equation 

H ^ cH + H c 

where = stands for the equality in distribution, < c < 1 and H c is a random variable inde- 
pendent of H. It is well-known that the law of these random variables is absolutely continuous, 
see e.g. [331 Example 27.8]. Moreover, Wolfe [36J showed that the density h of a positive self- 
decomposable random variable is unimodal, i.e. there exists a E 5K + (the mode) such that h is 
increasing on ]0,a[ and decreasing on ]a, oo[. The Laplace exponent, <j> s , of a self-decomposable 
distribution concentrated on 9 £ l + is given by 

oo 



(j) s {u) = au+ I (1 - e- UT )^dr 



o 



where A; is a positive decreasing function. We refer to the monographs of Sato [33] and Steutel 
and van Harn [33] for an excellent account on these sets of probability measures. 



1.1.2. Levy and Lamperti. Let ¥ x (we write simply P for Pq) be the law of a spectrally negative 
Levy process £ := (£,t)t>o, starting at x E with (F t )t>o its natural filtration. This law is 
characterized by the characteristic exponent of £i, which we assumed to belong to CIC, i.e. being 
of the form (jl.2p . We deduce, from the above discussion and the strong law of large numbers, 
that lim^ +00 ^ = sgn(E[£i])oo a.s. and the process oscillates if E[£i] = 0. 
For any 7 > 0, we write for the law of the Levy process with characteristic exponent 

V> 7 (u) = ip(u + 7) - ^(7), u > 0. 
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The laws pM and P are connected via the following absolute continuity relationship, also known 
as Esscher transform, 

(1.4) dP&V* = ert*- x '*-1'W t dP x \ Ft , t>0,xeK. 

Lamperti [23] showed that there exists a one to one mapping between ¥ x and the law Q e x of 
a ^-self-similar Markov process X on (0, oo), i.e. a Feller process which enjoys the following 
ce-self-similarity property, for any c > 0, 

(1.5) ((x cat ) t >o,® ( S) = ({cx t ) t > ,q^ 



More precisely, Lamperti showed that X can be constructed from £ as follows 
(1-6) log(X t ) = U t ,t>0, 



where 



A t = inf{s > 0; S s := f e aiu du > t}. 
Jo 



We write E x (resp. E x ) the expectation operator associated with (resp. Q x = Q x )■ 

Moreover, for E[£i] < 0, it is plain that X has a a.s. finite lifetime which is kq = inf{s > 0; X s - = 
0, X s = 0}. However, under the additional condition < 6 < a, where we recall that ift(0) = 0, 
Rivero [32], showed that the minimal process (X,kq) admits a unique recurrent extension that 
hits and leaves continuously a.s. and which is a a-self-similar process on [0, oo). We simply 
write (X,Q X ) for the law of such a recurrent extension starting from x > 0. Furthermore, for 
E[£i] > 0, Bertoin and Yor [5j Proposition 1], showed that the family of probability measures 
(Qx)ie>o converges in the sense of finite dimensional distribution to a probability measure Qq+ as 
x — > 0+, see also Caballero and Chaumont [10] for conditions for the weak convergence. Thus, 
for any x > 0, (X,Q X ) is also spectrally negative, in the sense that it has no positive jumps. 
Moreover, for any x > 0, (X,<Q X ) is a Feller process on [0, oo) and we denote its semigroup 
(resp. its resolvent) by (Qt)t>o (resp. by U q ,q > 0), i.e. for any x,t > and v £ 2?([0, oo)), the 
space of bounded Borelian functions on [0, oo), 



Qtv(x) = E x [v(X t )] 

/•oo 

U q v{x) = / e- qt Q t v(x)dt. 



o 



We also introduce the semigroup and the resolvent of the minimal process (X, kq), for x > 0, 

Q° t v(x) = E x [v(X t ),t < «o] 

/•oo 

U«v(x) = / e- qt Q° t v(x)dt. 



The strong Markov property yields the following expression for the resolvent of the recurrent 
extension 

(1.7) U q v(x) = U^v{x) + E x [e- qKo ]U q v(0), x > 0. 

Moreover, we recall that the Lamperti mapping reads in terms of the characteristic operator 
as follows. 
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Proposition 1.1. Let f : £K + — ► 9\ be such that f(x),xf'(x) and x 2 f"(x) are continuous 
functions on , then f belongs to the domain, B(L), of the characteristic operator L of(X, Q) 
which is given, for x > 0, by 



(1.8) l/(x) = x- a dx 2 r(^)+ter(x) 

+ ((/(e-^x) - /(x)) + xf(x)r) v{dr)^j . 

In the case E[£i] < 0, / £ D(L) i/ and onZy if it satisfies the boundary condition 

(1.9) lim = 0. 

x^O X y ~ L 

Proof. The first part of the proposition follows from Lamperti |23[ Theorem 6.1]. We point out 
that, in the case E[£i] < 0, the characteristic operator of the minimal process and its recurrent 
extension coincide for x > 0. Indeed, from its Feller property, the semigroup (Qt,t > 0) 
corresponding to the recurrent extension leaves invariant Co([0, oo)), the space of continuous 
functions vanishing at infinity. Therefore, if we let (|1.8p stand for the characteristic operator of 
the process, then we can conclude, see Gikhman and Skorohod [151 P- 130, Theorem 1], that the 
domain D(L) of the strong infinitesimal generator, L, of the process (X, Q) consists of all the 
functions / £ Co([0, oo))nB(L) such that L/ £ Co([0, oo)). However, it is plain that the analytic 
form of the function L/(x) for x > and for any function /, such that f(x),xf'(x),x 2 f"(x) 
are continuous, does not depend on the method of extension and is given by the expression 
above. Let us now turn to the boundary condition in the case < 0. We first deal with 

the necessary condition. From the discussion above, it is clear that it is enough to characterize 
the boundary condition for the strong infinitesimal generator. To this end, we recall that U q is 
a Fellerian resolvent, see Rivero |32| Theorem 2]. Thus, let / £ Co([0, oo)) (1B(L) then there 
exists v £ Co([0, oo)) such that L/ — qf = v, i.e. for x > 

U q v{x) = f(x). 

Then, using the expression of the resolvent of the recurrent extension (jl.7p and the continuity 
of /, we get 

(1.10) f{x) - /(0) = U*v(x) - E x [l - e-« Ko ]U«v(0). 

Let us now consider vo £ C*o((0, oo)), the space of continuous function vanishing at and 
oo. Then, from e.g. (32J Lemma 1], we have U g vo(x) £ Co((0, oo)) D O(L). Moreover, from [321 
Theorem 2] , we know that a necessary condition for the existence of a unique recurrent extension 
which hits and leaves continuously a.s. is that both limits 



*0 X y x^0 x y 

exist for any vq £ Co((0,oo)). Moreover, in this case, the identity 
(1.11) lim U ° Vo{x) - = U«v (0) 

holds for any vo £ Co((0, oo)). Hence, for vq £ Co((0,oo)), the necessary condition (|1.9p follows 
by dividing both sides of (jl.lOp by x e , by taking the limit x — > and by invoking the uniqueness 
of the limit. The general case, i.e. for any v £ Co([0, oo)), follows from the fact that the 
Lebesgue measure of the set {t >0;Xt = 0} is, by construction, with probability 1, and from 
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Blumenthal [HJ Section 4]. The sufficient part is readily obtained from the uniqueness of the 
recurrent extension that hits and leaves continuously a.s. □ 

1.1.3. The family of power series. Let ip G CKL and for 7 > and a > 0, set 

a n (V> 7 ; a) = JJ ip 7 (ak) , a = 1, 



vfc=l 



where we recall that ^(u) = ip(u + 7) — ^(7), ii > 0. Then, we introduce the function I ( 
which admits the series representation 

00 

T a ^(z) = y^a w (y> 7 ;g)z", z£t 



n=0 



We simply write T a r p when 7 = 0. We gather some basic properties of this family of power 
series which will be useful for the sequel. 

Proposition 1.2. For any ip E CK, 7 E £, 9^(7) > and a > 0, is an entire function. 

Moreover, i/E[£i] > or i/E[£i] < and 6 < a, T a ^ is positive and increasing on [0, 00). 



Proof. Observe that 



«n+l 



1 



I a n I I + 7) - V'(t) I ' 

The analyticity of 1 a ,ifa follows from the fact that lim u _ s . 00 ipy(u) = +00. The positivity and the 
monotonicity is secured by observing that under the condition 6 < a in the case E[£i] < 0, we 
have ip{ot) > and ip is increasing on [max(#, 0), 00). □ 

Remark 1.3. Note that Rivero's condition, < 9 < a for E[£i] < 0, arises naturally in the 
previous proposition to ensure that the associated functions are positive and increasing. 



2. Main results 

Let i/j E CKL. Moreover, if E[£i] < 0, we assume that 8 < a, recalling that ip(0) = 0. Next, for 
o G !EH, we introduce the stopping times 

r a = inf{s > 0; £ s = a} and K e a = mf{s > 0; X s = e a } 

with the convention that inf{0} = 00. For any A > 0, we denote p = 0(A) where eft : [0, 00) — * 
[max(#, 0), 00) is the increasing and continuous inverse function of ip. 

Theorem 2.1. Let q > and < x < a. Then, we have 



(2.1) 

Moreover, for A > 0, 

(2.2) E, 
and 

(2.3) E a 



E x [e 



-qKa 



e - qKa -XA Ka][{Ka<Ko} 



e J1 {r a <+oo} 



^a,i/.(g^°) 

T a ^{qa a )' 
'x\p T a ,i!> p 



(qx c 



In 



p(x-a)j^Pp_ 



2a,^ (qe™) ' 



Remark 2.2. (1) Note, by letting q -> in ([23]) . that [/s a < +oo] = 1 for any < x < a. 
Thus, the points above the starting point are recurrent states for (X, Q). Moreover, for 
> and x < a, we recall that P x [r a < +oo] = 1 and Q e ^ [n e a < kq] = 1. Thus, 
under such a condition, the indicator functions in (12.21) and in (12.31) can be omitted. 



(2) Consider the case ip(u) = \u 2 + ju, 7 > —1, i.e. (X, Q) is a Bessel process of index 7. 
Then, the expression of the Laplace transform of K a , which is well-known to be expressed 
in terms of the modified Bessel functions of the first kind, dates back to Ciesielski and 
Taylor [11] and Kent [20] . We refer to Section 14.11 for more detailed computations related 
to the modified Bessel functions. 

We proceed by characterizing, through its Laplace transform, the law of the exponential func- 
tional, Sqo, of spectrally negative Levy processes satisfying Rivero's condition. We emphasize 
that it is a companion result of Bertoin and Yor [51 Proposition 2] who computed the negative 
entire moments of the exponential functional of spectrally positive Levy processes when they 
drift to —00. 

Theorem 2.3. Assume E[£i] < and < 9 < a. Then, there exists a positive constant Cg 
such that 

T a ,ip {x a ) ~ Cgx e l a ^ g (x a ) as x —* 00. 

Moreover, if we assume that there exists f3 E [0, 1] such that lim^^oo ip(u)/u 1+/3 = lp then we 
have 

where M 7 = 0.577. . . stands for Euler-Mascheroni constant and a = ^ < 1. 
Next, introduce the function 

N"a,il>,o( x ) = ( x ) ~ C e x«l a ^ e (x) , x > 0. 

Ma,tj},9 is analytical on the right half plane and decreasing on 9 £ l + . Finally, the positive random 
variable has the following Laplace transform 

(2.4) E[e-« s °°] = M a ^ e (q)- 

Remark 2.4. Note that the random variable T,^ is solution to the random equation, for any 
a > 0, 

Soo = S r _ a +e 

where is an independent copy of Sqq. Indeed, together, the strong Markov property, the 
stationarity and independency of the increments of the Levy process £ entails that, for any a £ 9\, 
the shifted process (£t+ Ta — S,r a )t>o is distributed as (£t)t>o and is independent of (£t,t < r a ). 
Finally we get the equation by noting that, since E[£i] < 0, we have for any a > 0, P(r_ a < 
00) = 1. 

Bertoin and Yor [5] determined, in terms of their positive entire moments, the entrance law 
of spectrally negative self-similar positive Markov processes when E[£i] > 0. In the sequel, 
we characterize the entrance law of the dual process of (X, Q) when —00 < E[£i] < 0, i.e. of 
spectrally positive self-similar positive Markov processes when the underlying Levy process, in 
the Lamperti mapping, has a finite positive mean. To this end, let (X, Q) be the self-similar 
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process associated to the Levy process (£, P), the dual of (£,P) with respect to the Lebesgue 
measure. We recall that for — oo < E[£i] < 0, Bertoin and Yor [U Lemma 2] showed that, for 
x > 0, (X,Q X ) is weak duality, with respect to the reference measure m(dy) = ay a ~ 1 dy, with 
the minimal process (X, kq). In the same vein, Rivero |32^ Lemma 7] proved that in the case 
E[£i] < and 9 < a, (X, Q) is in weak duality, with respect to the measure m e (dy) = y a ~ e ~ 1 dy, 
with (X, Q^), the unique recurrent extension which hits and leaves continuously a.s., of the 

self-similar process associated, via the Lamperti's mapping, to (£ e ,P), the dual of the #-Esscher 
transform of (£,P). Before stating the next result, we recall that an entrance law {rj s ; s > 0} for 
the semi-group Qt is a family of finite measures on the Borel sets of (0, oo) such that r] s Qt = r/ s+t 
for all strictly positive s and t and such that [1 — e~ K °] remains bounded as s approaches 0. 

Corollary 2.5. If —oo < E[£i] < 0, then (X, Q) admits an entrance law which is absolutely 
continuous with respect to the reference measure m{dy). Its Laplace transform with respect to 
the time variable is given, for y, q > 0, by 

n q {y) = -^^Na^eiqy^. 



Moreover, assume E[£i] < and < 6 < a. Then, (X, Q' ') admits an entrance law which is 
absolutely continuous with respect to the reference measure m e {dy). Its Laplace transform with 
respect to the time variable is given, for y,q > 0, by 

fifty) = ^c e ^,e{ q y a ). 

Finally, we show that some specific combinations of the functions I a ,ib define some mappings 
from the convex cone CK, into the convex cone of positive self-decomposable distributions or 
into the convex cone of positive infinitely divisible distributions. 

Theorem 2.6. Let q > 0. Then, the mapping 

1 



(2-5) , T , , 

is the Laplace transform of a positive self-decomposable distribution. 
The mappings 

, 9fi , ( , i ( 

(2.6) q i ^ exp I -q-j- — -— and q ^ - — -- exp -q— — — 

are the Laplace transforms of positive infinitely divisible distributions. 

Finally, for E[£i] > 0, A>0 ; 0<a<^4<oo and recall that p = 4>(X), the mapping 

(9 ?) x / ayT a4 , p (a)l a ^{A) 



-AJ T a ^ p (A)I a ^(a) 
is the Laplace transform of an infinitely divisible distribution on the positive line. 

Remark 2.7. The random variable which Laplace transform is given in (|2.T[) is characterized 
below in (|3.7p . Moreover, consider again the case ip(u) = \u 2 + 7«, 7 > 0, i.e. (X, Q) is a Bessel 
process of index 7. 



(1) Then, (|2.T|) corresponds to Hartman's result f)l . 1 j) . Moreover, by letting A —* oo and 
using the asymptotic behavior of the modified Bessel function of the first kind, see (|4.ip . 
we get the Laplace transform of the so-called Hartman- Watson law [T7], the density 
of which has been characterized by Yor [38J. Note also that this law is the mixture 
distribution in the representation of the Von-Mises distribution as a mixture of wrapped 
normal distributions, see |17j . 

(2) Moreover, by choosing 7 = 5, the mapping on the right hand side of (|2.6p corresponds to 
the Laplace transform of the Levy's stochastic area integral. Indeed, for Bt = (Bj, B 2 ), 
a Brownian motion on 9\ 2 , Levy [25j computed the Laplace transform of the process 
L t = J* B]dB 2 s - B 2 s dB],t > 0, for fixed u > 0, and a = (s/u, y/u) G Dv 2 , as follows 

tu 

E[e ttLu \B u = a] = — — — exp -(tu coth(tn) - 1), t € 9*. 
smh(tn) 

This has been generalized by Biane and Yor [7] to the Levy's stochastic area integral 
associated to some planar Gaussian Markov processes, in terms of the modified Bessel 
functions of the first kind of any index 7 > 0. 



3. Proofs 

3.1. Proof of Theorem 12.11 First, since the mapping x 1— ► I a ^(x a ) is analytic on the right- 
half plane, it is plain that I a ip G B(L). Observe also that, for any > 0, x@ G B(L) and 

(3.1) Lx? = x p - a i/j(P). 

Then, for any positive integer N, using the linearity of the operator L and (|3,ip . we get, for any 
x > 0, 

N N 
Lj2an(^a)q n x an = ^a n {^a)q n I,x an 

n=0 ra=l 

JV 

r a(n—l) 



^2a n (Tp;a)q n Tp(c 



(ocn)x 
n=l 
N-l 

= q^2a n (^a)q n x an . 

n=0 

The series being analytic on the right-half plane, then the right-hand side of the previous line 
convergences as N — > 00. Hence, we get by monotone convergence (the series has only positive 
terms) that 

Ll a ^(qx a ) = ql a ,ip(qx a ), x>0. 
Moreover, recalling that for E[£i] < 0, 9 < a, we derive, in this case, that 

lim x -6+l°z a ( x a) = Urn V a(n + l)a n+1 (V, a)x a(n+1) - e 

n=0 

= 0. 

Hence, for E[^i] < 0, satisfy the condition (jl.9p . Thus, it is an eigenfunction for the 

recurrent extension (X, Q). Next, applying Dynkin's formula [12] with the bounded stopping 
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time t A n a , we get, for any t > 0, 

E x [e-^ AK «h a ^(qX? AK j] = l a ^(qx a ). 

Since the mapping a ^ X a ^(qa) is increasing on [0, oo), we obtain, by dominated convergence 
and by using the fact that the process has no positive jumps, that, for any < x < a, 

The proof of (|2.ip is completed. Next, the Esscher transform (jl.4p combined with the Doob 
optional stopping theorem yields 

e i {T a <+oo} — e fi i [ e J1 {r a <+oo}J • 

The proof of Theorem 12.11 is completed by invoking the obvious identity (K e a, Q e z) = (E Ta ,f x ) 
on {n e a < kq} and by using (|2.ip . □ 

We end up this part by providing an interesting absolute continuity relationship between 
self-similar processes, which will be useful for the sequel. We recall from the discussion above 
that for E[£i] < 0, the self-similar process associated to (£,P) has a a.s. finite lifetime Kq. Then, 
by using the Lamperti's mapping (jl.6p and applying the chain rule, it follows that the Esscher 
transform (jl.4p . reads, for j,5 > and x > 0, as follows 

(3-2) d®% t = (f) 7 " e-(^ *j« v on {i < « o}) 

where 

ft 

A = / X7 a du. 



o 

Since for every -stopping time T, £(T) is a F t -stopping time, the absolute continuity rela- 
tionship (|3.2p holds for every -stopping time on i*A T+ U {T < kq}. We mention that, for 
E[£i] > 0, where the condition "on {t < kq}" can be omitted, the relationship (|3.2p was already 
established by Carmona et al. [HI Proposition 2.1], under the name of Girsanov power transfor- 
mation. Finally, if we set 7 = 9 and 5 = in the case E[£i] < 0, then (|3.2p simplifies to the 
following Doob's /i-transform, for x > 0, 

(6) _(X t 



(3.3) dQ^ = ( ) dQ x , FAt , on {t < k } 



x 



3.2. Proof of Theorem 12.31 We assume that E[£i] < and 9 < a. Then, the identity 
(x a £oo,P) ( = } (ko.Q.) yields 

E x [ e - gK0 ] = E [e- qxa ^°°] . 

Hence, the mapping x \— > E x [e _<?K °] is decreasing on [0, 00), and by dominated convergence we 
have, see also Vuolle-Apiala [35] . 

lim E x \e- gK0 ] = lim E \e- qx "^] 

(3.4) = 

10 



and 



lim E x [e 



x->0 



We now compute the Laplace transform of kq. From the Doob's /i-transform (|3.3h . we deduce 
that, for any x, a such that < x < a, 

•a J l a ^ g (qa a )' 

Then, the strong Markov property and the absence of positive jumps yield 

1 



E, [e-^I W < Ko} ] 



E x [e~^l {K0<Ka} ] 



E [e~i Ka ] 



-qKa 



l a ,i(, (qx a ) - T a $ (qa a ) (~J 



E x [e~^I {Ka<Ko} ]) 
/x\6 la^g {qx a ) 



a J 1, 



{qa a ) 



and 



E* [e 



E x [e-« K °I {K0<Ka} ] +E X [e-^J {Ka<Ko} ] E a [e 



a \ / x\® 1 a ^ g {qx 



+ ~ 



E a [e 



a J Za^ g (qa a )' 
I a ,i> (qx a ) - x°l a ^ e (qx a ) 



(l^ (qa a )-E a [e-*»°]). 



To derive an expression of the sought quantity, we first differentiate with respect to a the previous 
equation, the function involved being smooth, to get the Riccati equation 







d 



—l a ^ (qa a ) ) + (E a [e~^] - (qa a )) 







■,-8 







dal a rf e (qa a ) 

where we have used the fact that x e l a ^ g (qx a ) > for any x > 0. By solving this Riccati 
equation, we observe that the solution has the following form 

(3.5) E x [e~^} = AI Q ^(qx a )-Cx e l a ^ e (qx a ) 

for some constants A, C. It is immediate that A = 1 since Eo [e™" 5 " 10 ] = 1. Then, the self- 
similarity property yields C = q&c, for some real constant c. Furthermore, (|3.4p ensures the 
existence of a constant Cg > such that 

(x a ) ~ C e x d l a ^ g (x a ) as X — ► oo. 

Moreover, observe, from (|3.5p . that 

E x [1 - e- K °l 



lim 



Cg. 



Next, we recall the following identities, see |32l (11) and Remarks 1. p. 489], 

E x [1 - e -"°] £(1-|) L^-il 

hm = TTT^r^ 

x^o x e aip 1 (0) 

li 



and 



E 



E< fl > 



-at. 



ds 



2—1 



Since (£,P' ') has a positive mean, the proof of the Theorem 
sition 2.3 in Maulik and Zwart [26] and by choosing c = Cq. 



is completed by using Propo- 

□ 



3.3. Proof of Corollary 12.51 We characterize the entrance laws of (X,Q) and (X, 
this end, we state the following easy result. 



. To 



Lemma 3.1. Let us assume that — oo < E[£i] < 0. Then, the entrance law of (X,Q) admits 
a density r]t(y) with respect to the reference measure m(dy), y > 0. Moreover, the Laplace 
transform in time of r] t (y), denoted by n g (y), q > 0, is characterized by the identity 



n q {y) 



E [c 



Proof. The claims follow readily from Bertoin and Yor [5, p. 396]. Indeed, they characterize the 
g-potential of the entrance law of a self-similar process associated, via the Lamperti's mapping, 
to a Levy process with a positive and finite mean, as follows. For any measurable function 
/ : 9\ + -> we have 



n q f 



|E[£ 



1 f°° 

ml mE[l 



m(dy) 



where we have used the identity E[£i 



-E[6 



□ 



The first part of the Corollary follows from (|2.4p . For the second statement, we use 
Proposition 3] where the g-potential, n q , of the entrance law is given, for a bounded continuous 
function /, by 

1 



/(y)E [e 



,a-l- 



y 



1 dy 



f%C 6 j 

The proof of the Corollary is completed by means of the identity 

Note the following result regarding the resolvent of the minimal process (X, kq). 

Lemma 3.2. For E[£i] < 0, the resolvent, Uq, of the minimal process (X,ko) admits a density 
with respect to the reference measure m{dy), which is jointly continuous and bounded. Thus, the 
processes {X,kq) and (X,Q) are in classical duality. 

Proof. We recall that in the case E[£i] < 0, the process (£,P) is necessarily of unbounded 
variation since we have excluded the case of negative subordinators. Thus, each point of the 
real line is regular for itself, see O VII, Corollary 5]. Moreover, since the g-capacity of {0}, 
which is (f>'(q), is positive for any q > 0, see O VII. 5. 2], and the resolvent of (£, P) is absolutely 
continuous, we deduce that points are not polar for (£,P), i.e. P x (Ty < +oo) > for any 
x,y £ 9\. It is not difficult to see that these two properties are left invariant by time change 
with a continuous additive functional, see Bally and Stoica [21 Proposition 4.1]. The assertions 
follow from [21 Proposition 3.1]. □ 
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3.4. Proof of Theorem 12.61 The claim (|2.5p is contained in the following Lemma. 



Lemma 3.3. The process (K a )a>o is under Q + an a- self- similar additive process, i.e. a process 
with independent increments which enjoy the scaling property (|1.5p . Hence k\ is under Q + a 
positive self-decomposable random variable. 

Proof. The first assertion follows from the absence of positive jumps, the strong Markov prop- 
erty and the self-similarity of (X, Q). The last statement is a straightforward consequence of 
the property of the law of additive processes, see Sato [33j Chap. 3.16] . □ 



Moreover, it is well know, see Wolfe [37] and also Jeanblanc et al. [19] for related results, 
that if the mapping q i— ► f(q) is the Laplace transform of a positive self-decomposable ran- 
dom variable then there exists a unique, in distribution, (increasing) Levy process L such that 
E [log(l + Li)] < +oo and its Laplace exponent, <f>L, is given by 

* i \ ^ Kq) > n 

Thus, from Theorem 12.31 and Lemma \3.3\ we deduce the statement (j2.6|) . after recalling that the 
infinite divisibility property is stable under convolution. 

Before stating the second Lemma, we assume E[£i] > and we introduce some notation. 
Let (Pt)t>o be the semigroup of the Levy process £. We denote by (Pp)t>o the subordinate 
semigroup of (Pt)t>0 by the continuous decreasing multiplicative functional (e~ q f )t>o- That is 
for / £ B(D\), we have, for any t > 0, 

Pff{x) = E x [e~^f{it)] , x £ JR. 

Next, by choosing A = in f)2.3f) . we deduce that the function x \— > (qe ax ) is excessive for 
the semigroup (Pp)t>o- Moreover, it is plain that, for any x £ SR, < T a ^ (qe ax ) < oo. Thus, 
one can define a new real-valued (sub)-Markov process with semigroup (resp. law) denoted by 
{Pf)t>Q (resp. P x ), as Doob's /i-transform of (P^)t>o, as follows, for any / £ B(d\) and t > 0, 

(3-6) Pff{x) = j \ qeax) P? U 1 ^ (qe a -)) (x), x £ JR. 

We are now ready to state the following result which characterizes the random variable associated 
to the Laplace transform (|2.7p . 



Lemma 3.4. For any < x < a, A>0 and recalling that p = (f>(\) , we have 



(3.7) E; 



i 



e -\T a 



e' 



Proof. We deduce from (|3.6|) that the following absolute continuity relationship 

, s t e- q ^T aiJ) {qe a ^) 
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holds for any t > and x £ 9\. It is plain that this relationship remains valid on F T + n{T + < 00} 
for any Foe-stopping time T. Recalling that [r a < 00] = 1 for E[£i] > 0, we get that 



e -\Ta 



e ^AT a -gS T 



1<x,il> (qe ax ) 

l Q ,ip p (qe aa )l a ^(qe a *) 

where the last line follows from (12.31). □ 



Finally, it is plain from the absolute continuity (13.81) that the process under F 1 is also spectrally 
negative. Then, from the strong Markov and the absence of positive jumps, we have for any 
x < c < a, 

(r^) { ^ (r^) + (r a ,lF c ) 

where the random variable on the right-hand side are independent. Hence, (r a ,P^) is infinitely 
divisible. The proof of Theorem 12.61 is then completed. □ 

4. Some illustrative examples 

We end up by investigating some well-known and new examples in more details. 

4.1. The modified Bessel functions. We consider £ to be a Brownian motion with drift 
7 G 9^, i.e. ip(u) = \u 2 + 7« and we set a = 2. In the case 7 < 0, we have 6 = 27 and therefore 
we assume 7 > — 1. Its associated self-similar process is well known to be a Bessel process of 
index 7. In the sequel, we simply indicate the connections between the power series and 
the modified Bessel functions since the results of this paper are well known and can be found 
for instance in Hartman [16] . Pitman and Yor [30] and Yor [39]. We have 



k=l 



on T(n + 7 + l) 

r( 7 + 1) 



Thus, we get 



where 



X 2 ,^(s) = (x/2)-^ 2 r( 7 + 1)I 7 (V2x~ 
(x/2)^ +2n 



I 7 (x) = £ 



n=0 



n!r(7 + n + 1) 



stands for the modified Bessel function of index 7, see e.g. |24t Chap. 5]. The asymptotic 
behavior of this function is well known to 

(4.1) I-y( x ) ~ —j= as x ~* °°- 
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Thus, we obtain that 

AA 2iV)27 (x) = (x/2)^ 2 r(- 7 + 1) (l_ 7 (yfa) - I 7 (V2i 

= ( x /2)V2_1_k 7 (^ 

where 2K 7 (x) = T(l — 7)r(7)(I_ 7 (x) — I 7 (x)) is the MacDonald function of index 7. 

4.2. Some generalizations of the Mittag-Leffier function. In [29], the author introduced 
a new parametric family of one-sided Levy processes which are characterized by the following 
Laplace exponent, for any 1 < £ < 2, (3 > and 7 > 1 — g, 

(4.2) ^((3u + 7) - r/>( 7 ) = - + 7 - 1), - (7 - 1),) 

where (A;)^ = r ^j^ stands for the Pochhammer symbol. Its characteristic triplet is a = 0, 

"(«v) = 757^ — v: ^ — y < °' 

/3L(2 -£»)(!_ e v)e+i 

and 

6 7 = ^(7)^(7- i + e )-*(7-i)) 

where ^(A) = Y(Xy * s ^ ne digamma function. In particular, if 70 denotes the zero of the function 
7 6 7 , then for 7 > 7o 6 (1 - g, 0), E[£ x ] > 0. 

4.2.1. The case 7 = 0. P~2j) reduces to ^(u) = ±(it - l) e . Observe that (9 = 1 and t//(1) = 
Moreover, setting a = g, we get 

.1 T(g(n +1)-1) 

a n (g,0;g) = — , a = 1. 

T(g- 1) 

The series, in this case, can be written as follows 

X Qrte( x ) = T(g)S B)e {gx) 
Z e ,il>{x) = r(g- l)£ e:e -i(gx) 

where 

^r(fm + /3) 

stands for the Mittag-Leffier function of parameter g,/3>0. The function £ e ,o(z) was defined 
and studied by Mittag-Leffier [28J. It is a direct generalization of the exponential function. The 
generalization £ e> p{z) was given by Agarwal pQ following the work of Humbert |18j . A detailed 
account of these functions is available from the monograph of Erdelyi et al. [13] . Next, we recall 
the following Mellin-Laplace transform of the generalized Mittag-Leffier function, for A, (3 > 0, 
we have 

°° e-^ x xP-% tP (x<>)dx ~ ( A + 1 )^ 



(A + l)o-l 
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We deduce by invoking a Tauberian theorem, see Feller |14l XIII. 5], the following asymptotic 
behavior 

£ Q ,p{x Q ) ~ -e x x 1 ~^l(x g ) asx^oo, 
Q 

i 



with I a slowly varying function at infinity. Thus, C\ = — E 



-l 

ye 



grig-, 1 ) 3t ,h 
r(i-i) and 



Finally, we can state the following properties of the Mittag-Leffler function. 
Corollary 4.1. Let 1 < g < 2. The mappings 

i 

<?^7 — and Q ^ £ q,q-i(<i) ~ tZqAv) 
ZqAI) q - 1 

are respectively the Laplace transform of a positive self-decomposable distribution and a com- 
pletely monotone function. The mapping 

*-* eX P c ( \ 
is the Laplace transform of a positive infinitely divisible distribution. 
4.2.2. The general case. In this case, 

n 

an(P, 7; er 1 = II (O + 7 - - (7 - 1)„) , ao = i. 

k=l 

and write £ et p a (x) = S^Lo a n(A 7! £?)x n - We point out that this function is closely related to 
the power series introduced by Kilbas and Saigo [22] which have coefficients of the following 
form 



a n (P,i;g) 1 = Y[{g{f3k + <y) + a = 1. 



k=l 



From Theorem 12. 6| we deduce the following properties. 

Corollary 4.2. Let l<g<2, (3>0 and 7 > 1 — g. Then, the mapping 

1 



is the Laplace transform of a positive self-decomposable distribution. The mapping 



q 1— ► exp 



£0,0/3,7(0) 



is the Laplace transform of a positive infinitely divisible distribution. Finally, for < a < A and 
7 > 70; recalling that p = 4>(X), the mapping 



-A) £ Q>Q p a + p (A)£ e>Q p a ( a ) 
is the Laplace transform of a positive infinitely divisible distribution. 
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4.3. The power series associated to stable processes and a new generalization of the 
exponential function. Finally, we consider the Esscher transform of spectrally negative stable 
process, i.e. ip~/(u) = c g ((u + 7)^ — (7)^), 7 > 0, 1 < g < 2 and c Q > 0. Its characteristic triplet 
is a = 0, 

V y; T(2 - q) I y y ' y 

1 

and 6 = c e £7 e > 0. The inverse function of ip is <fi(u) = (c e u + 7 s ) e — 7 and 

n 

e 1 ^)" 1 = n ((a^ + 7) e -7 e )> a = 1. 
fe=i 

Such formulation motivates us to introduce a generalization of the factorial symbol, which we 
defined, for n £ 01 ,a £ <£, 9te(a) > and 7 G £, ^(7) > 0, by 

n 

(a, l) e ,n = n (( ka + f) c - and ( a > 7)<?,o = L 
fc=l 



Note the obvious identities 



(a,7)o,n = 
(a,7)i,« = 



(a,0) e>ri = a n n!(a,0) e _i jn 
(a, 7 ),,n = Wl,-) • 

V OL/ g,n 



Moreover, we have 



We write simply 



Observe that 



k=0 



Wv)-E (fl)7W - 



n=0 



\iml a i e (c e x) = e<*. 
eJ.i 



Corollary 4.3. Xei 1<£<2, ct>0 and 7 > 0. T/ien, i/ie mapping 

1 

a 1 — ► 

-^Q, 7, 

is i/ie Laplace transform of a positive self- decomposable distribution. The mapping 

( Q^JlV 
q i-> exp — — 

is the Laplace transform of a positive infinitely divisible distribution. Finally, for < a < A and 
writing p = 4>(X), the mapping 



A 



a\ ^o, M (a)^,o,g(^) 

Za,p, e ( A )Za,0, e ( a ) 
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is the Laplace transform of a positive infinitely divisible distribution. 
References 

1. R. P. Agarwal, A propos d'une note de M. Pierre Humbert, C. R. Math. Acad. Sci. Paris 236 (1953), 2031- 
2032. 

2. V. Bally and L. Stoica, A class of Markov processes which admit a local time, Ann. Probab. 15 (1987), no. 1, 
241-262. 

3. J. Bertoin, Levy processes, Cambridge University Press, Cambridge, 1996. 

4. J. Bertoin and M. Yor, The entrance laws of self-similar Markov processes and exponential functionals of Levy 
processes, Potential Anal. 17 (2002), no. 4, 389-400. 

5. , On the entire moments of self-similar Markov processes and exponential functionals of Levy processes, 

Ann. Fac. Sci. Toulouse Math. 11 (2002), no. 1, 19-32. 

6. , Exponential functionals of Levy processes, Probab. Surv. 2 (2005), 191-212. 

7. Ph. Biane and M. Yor, Variations sur une formule de Paul Levy, Ann. Inst. H. Poincare Probab. Statist. 23 
(1987), no. 2, 359-377. 

8. R.M. Blumenthal, On construction of Markov processes, Z. Wahrsch. Verw. Gebiete. 63 (1983), 433-444. 

9. Ph. Carmona, F. Petit, and M. Yor, Sur les fonctionnelles exponentielles de certains processus de Levy, Stoch. 
Stoch. Rep. 47 (1994), no. 1-2, 71-101, (english version in [§9] p. 139-171). 

10. M.E. Caballero, and L. Chaumont, Weak convergence of positive self-similar Markov processes and overshoots 
of Levy processes, Ann. Probab., 34 (2006), no. 3, 1012-1034. 

11. Z. Ciesielski and S.J. Taylor. First passage times and sojourn times for Brownian motion in space and the 
exact Hausdorff measure of the sample path. Trans. Amer. Math. Soc. 103 (1962), 434-450. 

12. E.B. Dynkin, Markov processes. Vols. I, II, vol. 122, Academic Press Inc., Publishers, New York, 1969. 

13. A. Erdelyi, W. Magnus, F. Oberhettinger, and F.G. Tricomi, Higher transcendental functions, vol. 3, McGraw- 
Hill, New York- Toronto-London, 1955. 

14. W.E. Feller, An introduction to probability theory and its applications, 2 nd ed., vol. 2, Wiley, New York, 1971. 

15. I.I. Gikhman and A.V. Skorokhod, The theory of stochastic processes, springer- verlag, new york-heidelberg- 
berlin ed., vol. II, 1975. 

16. P. Hartman, Completely monotone families of solutions ofn-th order linear differential equations and infinitely 
divisible distributions, Ann. Sc. Norm. Super. Pisa CI. Sci. IV-III (1976), 267-287. 

17. P. Hartman and G.S. Watson, "Normal" distribution functions on spheres and the modified Bessel functions, 
Ann. Prob. 2 (1974), 593-607. 

18. P. Humbert, Quelques resultats relatifs a la fonction de Mittag-Leffler, C. R. Math. Acad. Sci. Paris 236 
(1953), 1467-1468. 

19. M. Jeanblanc, J. Pitman and M. Yor, Self-similar processes with independent increments associated with Levy 
and Bessel processes., Stochastic Process. Appl. 100 (2002), 223-232. 

20. J. Kent, Some probabilistic properties of Bessel functions, Ann. Probab. 6 (1978), no. 5, 760-770. 

21. A. A. Kilbas and J.J. Trujillo, Differential equations of fractional orders: Methods, results and problems, Appl. 
Anal. 78 (2001), no. 1-2, 153-192. 

22. A. A. Kilbas and M. Saigo, On solution of integral equations of Abel-Volterra type, Differential Integral equa- 
tions 8 (1995), no. 5, 993-1011. 

23. J. Lamperti, semi-stable Markov processes, Z. Wahrsch. Verw. Geb. 22 (1972), 205-225. 

24. N.N. Lebedev, Special functions and their applications, Dover Publications, New York, 1972. 

25. P. Levy, Wiener's random functions, and other Laplacian random functions, Proc. Sec. Berkeley Symp. Math. 
Statist. Probab., Vol. II, Univ. California Press, Berkeley, 1950, pp. 171-186. 

26. K. Maulik and B. Zwart, Tail asymptotics for exponential functionals of Levy processes, Stochastic Process. 
Appl. 116 (2006), 156-177. 

27. P.A. Meyer, Processus a accroissements independants et positifs, Seminaire de probabilites de Strasbourg 3 
(1969), 175-189. 

28. G. Mittag-Leffler, Sur la nouvelle function E a (x), C. R. Math. Acad. Sci. Paris 137 (1903), 554-558. 

29. P. Patie, Law of the exponential functional of a family of one-sided Levy processes via self-similar continuous 
state branching processes with immigration and the Wright hypergeometric functions, Submitted (2007). 



18 



30. J. Pitman and M. Yor, Bessel processes and infinitely divisible laws, Stochastic Integrals (Proc. Sympos. 
Univ. Durham, Durham, 1980), Lecture Notes in Math. (D. Williams, ed.), vol. 851, Springer, Berlin, 1981, 
pp. 285-370. 

31. D. Revuz and M. Yor, Continuous martingales and brownian motion, 3 rd ed., vol. 293, Springer- Verlag, 
Berlin-Heidelberg, 1999. 

32. V. Rivero, Recurrent extensions of self-similar Markov processes and Cramer's condition, Bernoulli 11 (2005), 
no. 3, 471-509. 

33. K. Sato, Levy processes and infinitely divisible distributions, Cambridge University Press, Cambridge, 1999. 

34. F.W. Steutel and K. van Harn, Infinite divisibility of probability distributions on the real line, Monographs 
and Textbooks in Pure and Applied Mathematics, vol. 259, Marcel Dekker Inc., New York, 2004. 

35. J. Vuolle-Apiala, ltd excursion theory for self-similar markov processes, Ann. Probab., Vol. 22, No. 2 (1994), 
546-565. 

36. S.J. Wolfe, On the unimodality of L functions, Ann. Math. Stat. 42 (1971), 912-918. 

37. S.J. Wolfe, On a continuous analogue of the stochastic difference equation X n — pX n -i + B n , Stochastic 
Process. Appl., 12 (1982), no. 3 , 301-312. 

38. M. Yor, Loi de I'indice du lacet brownien et distribution de Hartman-Watson, Z. Wahrsch. Verw. Geb. 53 
(1980), 71-95. 

39. , Exponential Junctionals of Brownian motion and related processes, Springer Finance, Berlin, 2001. 

Institute of Mathematical Statistics and Actuarial Science, University of Bern, Alpeneg- 
gstrasse, 22, CH-3012 Bern, Switzerland 
E-mail address: patieOstat .unibe . ch 



19 



